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ABSTRACT 

Gravitational lensing of an optical QSO by a spiral galaxy is often counteracted 
by dust obscuration, since the line-of-sight to the QSO passes close to the center of 
the galactic disk. The dust in the lens is likely to be correlated with neutral hydrogen, 
which in turn should leave a Lya absorption signature on the QSO spectrum. We use 
the estimated dust-to-gas ratio of the Milky-Way galaxy as a mean and allow a spread 
in its values to calculate the effects of dust on lensing by low redshift spiral galaxies. 
Using a no-evolution model for spirals at z <^ 0.8 we find, in A = cosmologies, 
that the magnification bias due to lensing is stronger than dust obscuration for QSO 
samples with a magnitude limit B ^ 16. The density parameter of neutral hydrogen, 
fiHi, is overestimated in such samples and is underestimated for fainter QSOs. 



Subject headings: cosmology: theory - gravitational lensing — quasars: absorption 
lines 
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1. Introduction 

Gravitational lensing by spiral galaxies results in multiple images with sub-arcsecond 
separation. A spiral galaxy at a cosmological distance would lens a background QSO if the 
line-of-sight to the QSO passed within a few kpc from the disk center. Under these circumstances, 
the neutral hydrogen in the disk would inevitably leave a strong Lya absorption signature on the 
quasar spectrum. A significant HI column density was indeed inferred for the spiral galaxy lens 
B0218+357 which has an Einstein ring radius ~ 0."3 (Carrilli et al. 1993). 

Strong Lya absorption is expected to result even at large impact parameters from the centers 
of galactic disks (~ 10 — 20 kpc), where the effect of lensing is weak. In fact, damped Lya absorption 
systems with neutral hydrogen column densities ;> 10 20 cm -2 are identified routinely in the spectra 
of high redshift QSOs, with an abundance consistent with that expected for the progenitors of 
present-day disk galaxies (Wolfe 1995). The association between damped Lya absorbers and the 
progenitors of present-day disk galaxies is supported by several observational facts. First, the 
abundance of heavy elements at low ionization stages in the absorbers and the velocity field traced 
by them is consistent with values expected for disk galaxies (Turnshek et al. 1989; Wolfe et al. 
1993; Lu et al. 1993; Pettini et al. 1994; Lu & Wolfe 1994). Recent Keck HIRES observations 
by Wolfe (1995; see also Lanzetta & Bowen 1992) indicate that the weak low-ionization metal 
absorption lines in these systems often show the highest column-density component at either the 
red or the blue edge of the velocity profile. Such an asymmetry is characteristic of absorption by 
a rotating gaseous disk (Wolfe 1995). Second, observations of redshifted 21-cm absorption and 
emission from damped Lya absorbers indicate disk-like structures of galactic dimensions (Briggs 
et al. 1989; Wolfe et al. 1992), and Faraday-rotation observations are consistent with the existence 
of micro-Gauss magnetic fields in these systems (Wolfe, Lanzetta, & Oren 1992; Welter, Perry, 
& Kronberg 1984; Perry, Watson, h Kronberg 1993). Finally, the association of the damped 
Lya absorbers with star-forming galaxies is directly confirmed by direct imaging (Steidel et al. 
1994, 1995) in most of the cases. However, whereas some systems are clearly spirals, others show 
irregular morphologies (Le Brun et al. 1996). For simplicity, we focus only on spirals. A more 
general treatment of this problem would involve many more free parameters and is beyond the 
scope of this paper. Neverthless, we stress that if dwarf galaxies with relatively shallow potential 
wells make a significant fraction of all damped absorbers, then our analysis overestimates the 
importance of lensing. 

The mass density of neutral hydrogen at z « 3, SIhi ~ 0.005, is comparable to the mass 
density in stars today, and appears to steadily decline at lower redshifts (Lanzetta, Wolfe, & 
Turnshek 1995; Wolfe 1995 and references therein). The simplest interpretation of this result is 
that the material in present-day disk galaxies had already assembled at a redshift z ~ 3, while the 
observed gradual decline of flni(z) since then is simply a result of the consumption of the gas by 
star formation. 

Based on these indications and the evidence that the velocity width of the associated metal 
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lines is indicative of rotational velocities ~ 200 km s _1 (Wolfe 1995), Bartelmann k Loeb (1996) 
have adopted the view that damped Lya absorption in QSO spectra results from neutral hydrogen 
disks in the present-day potential wells of spiral galaxies. Following this assumption, they showed 
that damped Lya systems may lens their background QSOs. An important consequence thereof 
is that the statistics of the absorbers themselves, as inferred from the number and distribution 
of absorption troughs detected, can be affected by the magnification bias. Because of the 
magnification of their flux due to lensing, the number of QSOs observed above a given magnitude 
threshold is increased, and consequently the observed fraction of QSOs which show damped Lya 
absorption is enhanced. 

An effect which tends to counteract the magnification by lensing is obscuration of the source 
by dust. By now, there is evidence that damped Lya systems contain measurable quantities of 
dust. A dust-to-gas ratio of order 10% of the Milky- Way value is inferred from the enhanced level 
of reddening found in a sample of QSOs that have damped Lya absorption (Fall, Pei, k McMahon 
1989; Pei, Fall, k Bechtold 1991), and from the gas-phase depletion of Cr relative to Zn (see 
Meyer, Welty, k York 1989; Meyer k Roth 1990; Pettini, Boksenberg, k Hunstead 1990; Meyer 
k York 1992; Pettini et al. 1994; Wolfe et al. 1994; Smith et al. 1996; although Lu, Sargent, k 
Barlow (1996) challenge this view based on the abundances of other metals). By dimming the light 
of the background sources, dust reduces the number of QSOs above a given magnitude threshold 
(Ostriker k Heisler 1984; Fall k Pei 1988, 1993), thus opposing the effect of the magnification 
bias due to lensing. It is thus a relevant issue to consider these two effects simultaneously and 
to estimate their relative importance. This is the purpose of the present paper. The detection of 
lensing by low-redshift spiral galaxies, such as B0218+357 at z = 0.68 (O'Dea et al. 1992; Patnaik 
et al. 1993), 2237+0305 at z = 0.04 (Huchra et al. 1995), and B1600+434 (Jaunsen k Hjorth 
1997) implies that dust obscuration does not mask lensing altogether. 

The efficiency of gravitational lensing is expected to strongly peak in the redshift interval 
0-2 ^ z ^ 0.8. Aside from the tendency of cosmological lenses to maximize their influence around 
z ~ 0.5 (Schneider, Falco, k Ehlers 1992; Bartelmann k Loeb 1996), evolutionary phenomena 
such as the hierarchical merging of halos and the HI consumption by star formation, would tend 
to reduce the significance of lensing per HI column density at high redshifts. We therefore focus 
in this paper only on spiral lenses at low-redshifts z <J 0.8. At these redshifts, there are various 
indications that the potential wells of spirals did not evolve significantly due to merging: 

1. Substantial merging or accretion would have resulted in asymmetric distortion of the 
isophotes of spiral galaxies. Direct observations of the level of lopsidedness in spiral galaxies 
set an upper limit on the current accretion rate of companions with a mass of order 10% the 
parent galaxy mass, which is 0.25 companions per Gyr (Zaritzky k Rix 1996). This implies 
that the fractional increase in the mass of galactic disks due to major mergers is smaller 
than 15% since z ~ 0.5. 

2. The thinness and coldness of the stellar distribution in galactic disks imply that these disks 
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did not encounter substantial merging since gas was converted into stars in them. Toth & 
Ostriker (1992) have obtained an upper limit of 4% on the fraction of the Galactic mass 
inside the solar radius that could have accreted within the last 5 billion years, based on the 
observed scale and Toomre Q-parameter of the Milky- Way disk. 

3. Considerable evolution in the potential wells of spiral galaxies since z ~ 1 would have 
resulted in a considerable scatter in their Tully-Fisher relation due to variations in their 
formation histories. Eisenstein & Loeb (1996) have used Monte Carlo realizations of halo 
formation histories in hierarchical models of structure formation to limit any substantial 
evolution in the potential wells of spirals since z ~ 1, based on the observed tightness in 
their Tully-Fisher relation. 

4. The statistics of multiply imaged QSOs implies that galactic lenses (either ellipticals or 
spirals) did not encounter significant merging below z ~ 1 (Mao & Kochanek 1994; Rix et 
al. 1994). 

5. Measurements of the rotation velocities and luminosities of spiral galaxies at intermediate 
redshifts (Vogt et al 1996) imply that they are kinematically similar to nearby spirals. 

6. Deep imaging and spectroscopy in the field of 3C 336 (Steidel et al. 1996) show a lack of 
strong evolution in the population of galactic metal-line absorbers up to a redshift z ~ 1. 

For redshifts below unity, the observed evolution in Qui is also very weak (cf. Fig. 8 in Wolfe et 
al. 1995). Based on all of this evidence, we adopt a model in which the gravitational potentials 
and the HI+dust content of the local population of spiral galaxies had not evolved significantly 
between a redshift of 0.6 ± 0.2 and the present time. 

The paper is organized as follows. In §2 we summarize our parameterization of dust at high 
redshift, focusing on its relation to neutral hydrogen in damped Lya systems. We then examine in 
§3 the simple situation of gravitational lensing and dust obscuration of a single point source by a 
single lens. We turn to statistical considerations in §4, where we detail our method to calculate the 
combined effects of gravitational lensing and dust obscuration of a population of background QSOs 
by an intervening population of damped Lya absorbers. Finally, §5 summarizes our conclusions. 



2. Dust Model 

Direct evidence that damped Lya systems contain dust is provided by the reddening of 
background QSOs and by the depletion of Cr relative to Zn. The amount of dust is usually 
expressed in terms of the dust-to-gas ratio k by 
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where tb is the optical depth due to dust in the rest frame B band and N21 is the column density 
of neutral hydrogen in units of 10 21 cnr 2 . 

The reddening of background QSOs can be estimated by comparing the spectral indices of 
quasars with and without damped Lya absorbers. Typical dust-to-gas ratios in the damped Lya 
systems can then be inferred from this reddening estimate if a shape for the extinction curve is 
assumed. 

Another way to infer k is provided by measurements of the gas-phase abundances of Zn, Cr, 
and other trace elements. The underlying assumption is that, apart from overall scalings in the 
metallicity and dust-to-gas ratio, the relative composition of the interstellar medium is the same 
in the damped Lya systems as it is in the Milky Way (i.e., the abundance patterns and grain 
composition are the same in both). Measuring the ratio Zn/H and the depletion of Cr relative 
to Zn, one can estimate the fraction of gaseous Cr that has condensed into dust grains. Both 
methods of measuring k yield a typical dust-to-gas ratio for damped Lya systems at z ~ 2 that is 
roughly 10% of the Milky Way value k ss 0.8. 

The dust-to-gas ratio evolves with redshift due to various astrophysical processes, such as 
conversion of gas into stars, metal enrichment and inflows or outflows of gas from galaxies. 
However, since we focus on low redshift absorbers, we assume no evolution in their dust-to-gas 
ratio, 

k(z) ~ jfe(0) . (2) 

The dust optical depth r(z) in an absorber at redshift z is related to the optical depth in the 
B band tb by 

*>-««(:£;). (3) 

where £(A) is the ratio of the extinction at wavelength A to that in the B band. Observations 
(Mathis 1990; Savage & Mathis 1979) indicate that the extinction curve declines with wavelength 
roughly as £(A) oc A . For wavelengths not too far from the B band, a good fit is provided by 

i(X) ~ 0.4- fim . (4) 
A 

In terms of the dust-to-gas ratio as defined in equation (jl]), the optical depth can be written 

as 

r(z) = k(z)N 21 ^-^j . (5) 

The intrinsic flux S' of a source is changed to S if it passes an intervening galaxy. Gravitational 
lensing magnifies the flux by a factor fj,, while obscuration by dust reduces the flux by a factor 
exp(-r) (Fall & Pei 1989). The net effect is thus 



S = fj, exp(— r) S' 



(6) 
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3. Combination of Lensing and Dust Obscuration 

We now examine the combined effect of dust obscuration and gravitational lensing in the 
simple situation of a point source which is lensed and obscured by a single intervening galaxy. For 
definiteness, we adopt a source redshift of z s = 3 and a lens redshift of z& = z^s = 0.5 where 
necessary. 

We describe the profile of the neutral hydrogen disk in the galaxy by the radius where the 
face-on column density is N', 

Rr(N') = fu(N') RJ l , (7) 

where I = L/L* is the galaxy luminosity, t ~ 0.4, and R* ~ 11 h^ 1 kpc is a characteristic radius of 
an L* galaxy. The data by Broeils & van Woerden (1994) imply an exponential neutral hydrogen 
profile, 

^H-^ Gxlffon-O - <8) 

For simplicity, we will in this section consider a lensing galaxy with luminosity L*. The 
hydrogen radius of equation (|7|) then simplifies to Rk(N') = fn(N') i?*, where N' is the face-on 
column density of HI. If the disk is inclined by an angle 9 with respect to the line-of-sight, the 
observed column density is N = N'/ cos9. 

The intercepted HI column density not only depends on the relative position of absorber and 
source, but also on the lensing properties of the absorber because the light rays are bent. We 
model the lens as a singular isothermal sphere (SIS). Recent work (Kochanek 1996) has shown 
that, in the case of elliptical galaxy lenses, the addition of a small core radius, consistent with 
the kinematics and the image configuration of the lens, does not change the lensing probability 
significantly. Since the light distribution in the central bulge of spiral galaxies resembles that in 
elliptical galaxies, we assume the SIS model to be a good approximation for our purposes. 

In the SIS model, the scaled distances of the source y and of the images x from the center of 
the lens are related by (e.g., Schneider et al. 1992) 

x={ V±1 if ^ 1 (9) 
1 y + 1 if y > 1 v ' 



The appropriate length scale for x is 



V c J H r s 

where rd jS ,ds are the angular diameter distances from the observer to the lens, to the source, and 
from the lens to the source, respectively, in units of the Hubble length cHq 1 . In the numerical 
examples throughout the paper, we will use a Hubble constant Hq = 100 h kms -1 Mpc -1 , with 
h = 0.5, and an galaxy velocity dispersion of a v * ~ 160 kins" 1 . The magnification of a point 
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source by an SIS lens is 



J 2/» if 9 < 1 M1 , 
" = 1 + 1/9 «»>1 ' (U> 



In order to compute the optical depth due to dust obscuration, we need to calculate the 
column density N' of the galactic HI disk where it is intercepted by the light rays from a source 
at position y. Since an inclined disk appears elliptical, the radius in the plane of the disk where 
the light rays pass the galaxy depends on the relative position angle between y and the inclined 
galactic disk. For an inclination angle 6, the radius averaged over position angle is 



1 r r i 2x 
(x) = — / dip (x cos f) 2 + (x sin ipl cos 9) 2 = -E(sin 2 #) (12) 

2vT Jo L i IT COS O 

with the elliptical function of the 2 nd kind E(x). 

The HI column density N' in the disk at the average radius (a;) is then found from 

R*MN') = (x) (o (13) 

which can be solved for N' using fn(N') from equation (g). The observed column density N as a 
function of the source position y is 

N(y) = 2Xl f s g Cm " 2 lQ 3/4 " Ma)K ° /(2fl,) , (14) 

with x(y) given by equation (§). The optical depth can then be easily computed from equation 
For the redshifts chosen, z s = 3 and = 0.5, for observations in the B-band, and for a cosmological 
model with Qq = 0.2 ,$7^ = 0, we find from equations (Q) and (|l0|) 

Co ~ 4 kpc , £ m 1.4 , R* = 22 kpc , (15) 

for h = 0.5. 

We now have all the ingredients to estimate the ratio between the observed and intrinsic 
fluxes, S and S", related by equation (||). The resulting flux ratio S/S' is shown in Figure [j], where 
5/5' is plotted as a function of source position y for different values of the inclination angle 8 of 
the disk. For comparison, the solid line shows S/S' without taking dust obscuration into account. 
Figure [l] shows that for each inclination angle there is a source position such that magnification 
dominates over dimming (i.e., S/S' > 1) for y < y* while extinction is stronger than magnification 
(i.e., S/S' < 1) for y > y*. As smaller y imply smaller x, we find that magnification is stronger 
than extinction for light rays passing closer to the center of the disk. 



4. Statistics of Dust Obscuration and Lensing 



We now examine the combined effect of lensing and dust obscuration on a population of 
distant QSOs by a population of intervening galaxies. Damped Lya systems offer one of the 
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best ways of tracing the evolution of neutral hydrogen in spiral galaxies at intermediate and 
high redshifts. We therefore focus on the simultaneous effects of lensing and dust obscuration 
on the properties of neutral hydrogen which can be inferred from the damped absorption lines. 
In particular, we investigate the effects of dust and lensing on the column density distribution 
function of HI, f(N), and on the estimates of the contribution of HI to the density of the Universe 

We model the column density profile of neutral hydrogen in the absorbers according to 
equations (0) and (||) and assume for simplicity that it is constant in time. The absorber 
population is described by a Schechter luminosity function, 

$(/) dl = I' exp(-Z) dl , (16) 

for which we adopt the parameters appropriate for spiral galaxies (Marzke et al. 1994), 

= 1.5 x KT 1 h 2 Mpc~ 3 , s = -0.81. (17) 



The source population of distant QSOs is parameterized in terms of the intrinsic number 
count function which can be well modeled by a broken power law (e.g., Boyle, Shanks, & Peterson 
1988; Hartwick k Schade 1990), 

where the flux Sq at the break point corresponds to an apparent QSO magnitude of B ~ 19.5 and 
A is a constant. A reasonable fit to observational data for the two branches of the power law is 
provided by the values ft = 0.64 and ft = 2.52 (Pei 1995). 

If neither lensing nor dust obscuration are effective, the probability for a QSO to exhibit 
damped Lya absorption with a column density > N in its spectrum is given by (Bahcall & Peebles 
1969) 

P LYa (z 1 , Z2,N) = r(l + s + 2t)n Rl (&(N)) AX( Zl ,z 2 ) , (19) 

where AX(zi, Zz) is the absorption distance surveyed per QSO, and (f^(N)} is the average of (|8|) 
over all inclination angles. 

If the absorption probability PL ya (zi, Z2, N) is measured, the column density distribution of 
absorbers f(N) is given by 



/(AO ' ' 



dP(z l7 z 2 ,N) 



dN 



H AX(z 1 ,z 2 ] 

This expression is independent of the QSO redshift, the QSO flux, and the redshift interval 
surveyed for damped Lya absorption. 



(20) 



On the other hand, if gravitational lensing is taken into account, the probability for damped 
Lya absorption, and thus the column density distribution, does depend on the flux threshold 
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S and on the redshifts involved because of the magnification bias (Bartelmann & Loeb 1996). 
Specifically, 



/GL+dustOlj z 2 , z s ,N;S) = — 



1 



H AX(z u z 2 ) 



dz f G L +dust (z,z s ,N;S) 



(21) 



with 



fGL+du S t(z, Z S ,N;S) 



Nqso{S) Jo 



dS' 



dp' GL (z,z s ,N;S/S') 



ON 



(22) 



where p' GL (z, z s , N; fj,) is the probability that a QSO at redshift z s shows damped Lya absorption 
with column density > N and is magnified by a factor > // by lenses in the redshift interval dz 
around z (cf. equation (2.25) in Bartelmann & Loeb 1996). 

Let us now consider the combined effect of lensing and dust obscuration. For an intrinsic 
source flux S', the observed flux is S = S' fj, exp(— r) with r(k, z) = /c^i^As/l + z) (cf. eqs. (|5|) 
and Consequently, the distribution of absorbers in equation (|22| ) becomes 



/GL+dust(^, z s ,S;N) 



Nqso{S) Jo 
A 

N Q so(S) 



dS' 



dp* (z,z s ,N;S/S'e-T) 



So 



dS' 



ON 

0p' GL (z,z s ,N;S/S'e~ T ) 



+ 



dS' 



So 



dN 

dp' GL (z,z s ,N;S/S'e- T 



~dS^ [b} 



ON 



S' 



(23) 



where we have used equation (|18|). By changing variables S = S'e T , equation (|23|) becomes 



fGL+dnst(z,z s ,S;N) 



A 



+ e 



iV Q so(5) 

/3 2 r 



-Pit 



S e~ T 



dS 



dp' GL (z,z s ,N;S/S 



dN 



dS 



Soe- 



8p' GL (z,z s ,N;S/S 



ON 



S 



-/3 2 -i 



(24) 



Hence, dust causes an apparent shift in the break flux of the QSO distribution function and damps 
the absorber distribution by a factor exp(— (3t) (a discussion of the latter effect for bright QSOs 
can be found in Fall & Pei 1993). 

Since the value of k varies considerably among different DLA systems, we allow for a spread 
in its values. Following Fall & Pei (1993), we assume a log-normal form for the probability 
distribution of k, 



g(k) 



1 



27T<7* k 



exp 



(In A; — In k*Y 
2^! 



(25) 
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where the parameter k* is related to the mean k of the dust-to-gas ratio by 



k = dk g(k) k = k* exp ( — 
J o V 2 



(26) 



with k given by equation (||) under the assmption that at z = the Milky Way value of k reflects 
the true mean of the spiral galaxy population. Following the study of Fall & Pei (1993; cf. their 
figure 8b), we adopt a constant value of = 1. With the introduction of a probability distribution 



for k, equation (24) reads 

fGL+dust(z,Z s ,N, S) 



A 



dkg(k, z) 



Nqso(S) Jo 

-Ar(M f Soe ~ T<k ' Z) , , dp' GL (z,z s ,N,S/S>) , 
Jo 



dS' 



dN 



+ 



-/3 2 T(k,z) 



Soe -T(k,z) 



, dp' GL (z,z s ,N,S/S') rt-fa-x 
<lh dN b 



(27) 



Equations (^|) and (^7j) can be used to get the observed distribution of absorbers in the redshift 
interval Z\ < z a b s < z%. The results from these equations are displayed in Figure |2| for the absorber 
redshift interval 0.4 < z a b s < 0.8. This figure shows the observed distribution of N f(N) for 
two choices of the limiting flux B and for three cosmological models, in comparison to the true 
distribution (solid line). 

Having calculated Nf(N), the inferred column density of neutral hydrogen in absorbers with 
column density of N% < N < N2, normalized by the present critical density, is given by 



Hn ft) r N2 
^m{N x ,N 2 ) = dNNf(N) 

c p c ,0 JNi 



where m is the mean molecular mass per proton and p Cj o is the present-day closure density, 



Pc,0 



8ttG 



(28) 



(29) 



Figure || displays the result for Ohi for the absorber redshift interval 0.4 < z a b s < 0.8 and for the 
same three cosmological models as in Figure ||. 

The enhancement in the abundance of absorbers (and hence in the inferred density parameter 
of neutral hydrogen) is a consequence of the magnification bias due to lensing, which increases 
the number of QSOs above the detection threshold. Extinction by dust results in an opposite 
tendency. In Figure 4a we show the combined effect of these opposite tendencies as a function 
of the QSO magnitude limit, for low-redshift absorbers with iV ;> 10 20 cm -2 . The QSO detection 
probability is compared to the case where neither of these two effects are taken into account. We 
find that for bright quasars with a limiting magnitude B ^ 17, the effect of lensing dominates. 
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In Figure 4b we compare the above probability to the situation where only lensing is taken into 
account (cf. Bartelmann & Loeb 1996). This figure illustrates that for low-redshift absorbers, the 
suppression of the lensing effect by dust is of several tens of percent. The level of suppression is 
higher for brighter QSOs because of the dependence of the extinction factor on the slope of the 
QSO number counts [cf. eq. (p7|)], which is steeper for brighter QSOs. Note that the suppression 
of the lensing probability by dust is weaker than the suppression of the lensing effect on Qhi 
(compare Fig. 5 of Bartelmann & Loeb (1996) to Fig. || in this paper). This follows from the 
fact that Qui is integrated with an extra factor of iV (cf. eq. ^8|), which tends to emphasize high 
column density systems more, and therefore enhance the suppression of lensing by dust at a given 
QSO magnitude limit. 

5. Conclusions 

Gravitational lensing by a spiral galaxy occurs when the line-of-sight to a background quasar 
passes within a few kpc from the center of the galactic disk. Since galactic disks are rich in neutral 
hydrogen, the quasar spectrum will likely be marked by a damped Lya absorption trough at the 
lens redshift. Indeed, the observed 21 cm absorption in the spiral galaxy lens B0218+357 implies 
an HI column density of 4 x 10 21 cm -2 (T s /100 K)/(//0.1), where T s is the spin temperature 
and / is the HI covering factor in this system^]. Bartelmann & Loeb (1996) have shown that the 
efficiency of searches for gravitational lensing with sub-arcsecond splitting can be enhanced by 
1-2 orders of magnitude if one focuses on the subset of bright quasars which show low-redshift 
(-£abs ^ 0.8) strong Lya absorption (N ^ 10 21 cm -2 ). The double- image signature of lensing could, 
in principle, be identified spectroscopically and without the need for high-resolution imaging; the 
absorption spectrum of the lens might show a generic double-step profile due to the superposition 
of the two absorption troughs of the different images. The different images pass the absorbing disk 
at different impact parameters and therefore have different spectral widths. 

Lensing of optical sources by spiral galaxies is expected to be counteracted by dust 
obscuration, since the central opacity of nearby disks is considerable (see, e.g. Byun 1993; Davies 
et al. 1993; Byun, Freeman, & Kylafis 1994; Jansen et al. 1994; Rix 1995). The distribution of 
dust in galactic disks is correlated with the distribution of neutral hydrogen, which in turn can be 
probed in the absorption spectra of QSOs. In this paper, we used the estimated dust-to-gas ratio 
of our galaxy and allowed a spread in its values to calculate the effects of dust on lensing by low 
redshift spiral galaxies. For absorbers with z a b s <J 0.8, we have found that the relative importance 
of the magnification bias due to lensing and of the obscuration due to dust depends on the QSO 
magnitude limit and on the choice of cosmological parameters (cf. Fig. 4a). Overall, the effect of 
dust suppresses the detection probability of lensing by <J 55% (cf. Fig. 4b). 



lr The high HI column density in this lens will be directly calibrated based on the absorption spectrum of the QSO 
in forthcoming HST observations (Falco, Loeb, & Bartelmann 1996). 
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The magnification bias and the deflection of fight rays due to iensing change the column 
density distribution of damped Lya absorbers, f(N). Afthough dust tends to suppress the 
distortion to f(N) derived by Bartelmann k Loeb (1996), the generic iensing peak for low-redshift 
absorbers with iV ~ 10 21 ' 5 cm~ 2 is still evident in Figure 2a. This peak results from the 
magnification bias, which raises faint QSOs above the detection threshold, thereby enhancing the 
probability for observing damped Lya absorbers. The sharp cutoff at yet higher column densities 
arises because the bending of light rays due to Iensing yields a minimum value for the impact 
parameter of the line-of-sight to the QSO through the galactic disk. For low-redshift absorbers 
and models with no cosmological constant we have found that the density parameter of neutral 
hydrogen f^Hi; obtained by integrating Nf(N) over all column densities, is overestimated for QSO 
samples with B <; 16 and is underestimated for fainter QSOs (cf. Fig. 3). 

Since the current set of known damped Lya systems includes only a small number of 
low-redshift absorbers (see, e.g. Table 3 in Wolfe et al. (1995)), it is difficult to compare the 
statistical results shown in figures 2-4 to existing data. We expect the significance of Iensing 
relative to dust to decline at z ^ 1, and to yield only a secondary correction to the inferred 
statistics of the current set of damped Lya absorbers (Wolfe et al. 1995; Storrie-Lombardi et al. 
1996a, b). The quantitative predictions of this paper will be tested in a straightforward manner 
with the advent of absorption data for larger statistical samples of QSOs with uniform selection 
criteria, such as the set of ~ 10 5 QSOs in the forthcoming Sloan digital sky survey (Gunn k 
Knapp 1993). 
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Fig. 1. — The ratio between the observed flux S and the intrinsic QSO flux (without the absorbing 
lens), S' , as a function of the dimensionless impact parameter y. The solid curve shows the ratio 
where only lensing is taken into account, while the other two curves consider the combined effect 
of lensing and dust for a face-on disk (dotted line) and a disk inclined by 70° (dashed line). 
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Fig. 2. — Column density distribution of damped Lya absorbers per logarithmic column density 
interval, Nf(N), for samples of QSOs with flux limit B = 17 [panel(a)] and B = 19 [panel(b)]. The 
solid curve in both panels displays Nf(N) without the influence of gravitational lensing and dust 
extinction. The other three curves take into account both lensing and dust for three cosmological 
models characterized respectively by the parameters ^o=0.2, f2^=0.8 (dotted line); Qq=0.2, O,\=0 
(dashed line); and £lo=l, ^a=0 (dashed-dotted line).w In both panels the QSO redshift is fixed at 
z s = 3 and the redshift interval for the absorbers is 0.4 < z a b s < 0.8 
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Fig. 3. — The ratio between the inferred and the true S7hi for three cosmological models 
characterized respectively by the parameters Oo=0.2, fi,\=0.8 (dotted line); f2o=0.2, f2A=0 (dashed 
line); and f2o=l, ^a=0 (dashed-dotted line). 
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Fig. 4. — Probability for a QSO to be detected with apparent magnitude brighter than B and show 
damped Lya absorption with column density N > 10 20 cm -2 when lensing and dust obscuration are 
taken into account. The probability is normalized by its value in the absence of an absorber [panel 
(a)], or with only lensing taken into account [panel (b)]. In both panels, z s = 3, 0.4 < z a b s < 0.8, 
and the three curves refer to the same models as in Figures || and ^. 



